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Abstract. We consider the relation between so called continuous localization models - i.e. 
non-linear stochastic Schrodinger evolutions - and the discrete GRW-model of wave function 
collapse. The former can be understood as scaling limit of the GRW process. The proof relies 
on a stochastic Trotter formula , which is of interest in its own right. Our Trotter formula also 
allows to complement results on existence theory of stochastic Schrodinger evolutions by Holevo 
and Mora/Rebolledo. 



1. Introduction 

The Ghirardi- Rimini- Weber theory (shortly GRW, see [9], [3] and, for a mathematically rigorous 
treatment, [35]) is one of the best known versions of stochastic localization models, aiming at a 
solution of the measurement problem in quantum mechanics (see e.g. [5]). In the GRW theory 
the Schrodinger equation is interrupted by the random collapse of the wave function at randomly 
chosen times. In this way, one obtains a description of microscopic quantum phenomena where 
macroscopic superpositions (i.e. Schrodinger cat states) are surpressed. The GRW-theory contains 
two new constants of nature: the collapse of the wave function ip occurs at random times with 
intensity p and the wave function collapses in space by the multiplication with a Gaussian hitting 
function with spread 1/y/a. The centers of the Gaussians are approximately |^| 2 -distributed. 

A related class of stochastic Schrodinger evolutions, where the collapse of the wave function 
happens continuously, appears in the theory of continuous quantum measurement and in the 
theory of open quantum systems (jU H2 HH [T5j [H]). The Schrodinger equation is replaced by a 
Hilbert-space valued stochastic differential equation. 

A known example is the Diosi equation ( [7J [HI [201 HI 12] ) 

^ ^ dtpt(x) = —iHrptdt + VXip t {x)x ■ d£ t - ^\x\ 2 tp t (x)dt 

ipo = <t>0, 

where (£ t ) is a three-dimensional Wiener process on a filtered probability space (CI, {?, ({ft)) Q) an d 
A > 0. The physical wave function process arises from ip by normalization V'/ll^'ll an d by a change 
of measure. The normalized wave function obeys then a non-linear equation. 

In the physics literature one often considers the equations for the statistical operators pt '■= 
E(|^t)(^t|). The von-Neumann equation is then replaced by so called Lindblad equations (|17j). 
which for the GRW-process reads ([3]) 

j t Pt{x,y) = - l -[H,p t ]{x,y) -fi(l -e^l^l 2 ) Pt ( x ,y) 

and (0) 

—p t (x,y) = -^[H,pt](x,y) - ~\x - y\ 2 p t (x,y) 

for the Diosi equation. We observe that for \x — y\ 2 << 1/a and p ■ a = 2A both equations are 
approximately the same: 



M (l- e -fl*-*l 2 ) « M (l-(l- j\x-y\ 2 ) 



\\x-y?- 
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In [7] , Diosi addresses therefore the question, whether equation can be obtained as a scaling 
limit from the GRW-process by increasing the spread of the hitting function - a — ¥ - and, at 
the same time, the collapse frequency fj, — > oo , such that 

„, u • a 

1.2 = A = constant. 

v ' 2 

We shall prove that this is indeed the case. A natural method for the proof is suggested by 
the structure of the GRW-process itself: the Schrodinger evolution and the collapse mechanism 
act alternately. This has obvious similarities with the Trotter formula, which is well known and 
often used in many areas of mathematical physics. In our situation we need a stochastic version 
of the Trotter formula, since both the collapse times and collapse centers are random. Stochastic 
versions of the Trotter formula have been considered since long ([16), but the underlying state 
space has been assumed to be locally compact, which excludes infinite-dimensional Hilbert spaces. 
More recently, a version of the Trotter formula for a class of stochastic Schrodinger evolutions 
with deterministic jump times and bounded collapse operators has been established ([10). 

In the following section we introduce the GRW-process and the Diosi process. We observe that 
the collapse mechanism of GRW can be rephrased in terms of of increments of the Wiener process 
which appears in the Diosi process. This then allows to rephrase the convergence in terms of a 
Trotter product formula. In section 3 we first prove the Trotter formula. 



2. The scaling limit of GRW 
In the GRW-process the Schrodinger evolution is interrupted at the random times 

(2.D T - : =£f< 

fe=l ^ 

where X\, X 2 ■ ■ ■ are independent exponentially distributed random variables with EA& = 1. At 
the times T„. M the wave function gets multiplied by a random Gaussian function with spread -^j, 

where the centers of the Gaussians are random variables Y = (Yi , Y2, . . . ) the distribution of which 
is given by the wave function in the following way: 

We collect the randomness in the sample space fi := R N x (R 3 ) N with coordinate projections 
X = (Xi,X 2} . . . ) on the first and Y — (Yi, Y 2 , . . . ) on the second factor. We define P(r £ •) to be 
a countable product of exponential distributions, P(Y € -|r) will be specified via the wavefunction. 
Starting with (j>o £ L 2 (M 3 ,C), we define recursively 

(Yr, . • . , Y n ,X, x) := (-) 3 e"* l*- y »lV**» ff £r B _ 1 (Y T ^, . . . , Y Tn _^ , X, x) 

(2 ' 2) / 



and 
(2.3) 



\AY T ^ e A I Y Tlft ,..., Yr^ , X) := J || (Y Tl „ , . . . , Y Tn _^ , A, y 
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' ')\\L 2 (R3) 



dy 



J A \f^J e- alx - yl2 \(e-i x " H cl )Tn _ 1 )(Y,X,x)\ 2 dxdy, 



which reads explicitely as 
(2.4) 



'<*A Y Ti,» e4)= 11^ (j/i,A,-)||i2 (R 3 ) dj/i 



/ e -«l-f 1 | 2 |( e -*^ ff ^ )(x)| 2 dxdt/i, 
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: / P a ^(^, (1 e A 2 1 Y Tl ^X)dF a ^(- | X) 

• y { y T 1 . u eA 1 } 



{y Ti ^a 1 }Ja 2 V v J WnJ\ 
2 1 ' ^ [ e-^-y^\{e-i x > H ip Tl J( yi ,X,x)\ 2 dxdy 2 1 x J ^Jy^X^fdy 



y JTl j yi ,x,-w" yTl - 

a \*-y*\ 2 \( e -i X2H e-^-- yi \ 2 e-i XlH (f> )(x)\ 2 dxdy2dy 1 



i 



Ai J A 

f f [ " ,r 

J A! JA 2 7r 

and, inductively, 

P a , fl (Y Tliii eA 1 ,...,Y Tn ^eA n \X) 

(2 ' 5) = / J^/ Ke-fl-^^e-t^.-.e-fl-^^e-^^WI^d^,...,^). 

JA 1 x---xA„ V I" 7 

The sequence of wavefunctions is extended continuously according to 
(2.6) $■„,„+. e ~ isH ^T n>ll for < s < X n+1 

between the collapse times. By this and the preceding equation, the process {<j>t)t>o on the 
probability space (fi, P) is specified. 

In the Diosi process the wave function obeys the stochastic Schrodinger equation 

(2-7) dMx) = -iHiptdt + \f\^ t {x)x • d& - ^\x\ 2 Mx)dt, 

where H is the Hamiltonian which for the purpose of this paper can be thought of as the 
Schrodinger Hamiltonian H = — + V. The noise (£t) is a standard <i-dimensional Wiener 
process on a filtered probability space (fi, (3t), Q) and A a denotes positive intensity parameter. 
The "physical" collapse process is then given by 

<h ■= ^ 

\\lpt\\L2(M<i) ' 

weighted by a new measure defined as 

(2.8) F X (A) := E q (xa|!^|| 2 ) for A e &. 

The martingale property of \\ipt\\ 2 (see [IS]) ensures that this is indeed a proper definition of a 
measure on ft. 

We start the comparison of the GRW and the Diosi processes by first neglecting the Schrodinger 
evolution (i.e. we set H = 0) and the randomness of the collapse times (i.e. Xk = 1). Then, from 
(I2~2l and {23]), 

(2.9) (x) = (-Y 4 e -^-^i 2 +---+i*-^i 2 Vo(z), 



(2.10) 

and 
(2.11) 



||^ 



'AiX---xA„ 

The solution of the collapse part of (|2.7I) is 

(2.12) ^(x) = e^'^-^^'VoW, 
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and can be rewritten as 
(2.13) 

. , > EJUitVXai-Ki-CAri)-^!*! 3 ] , , , -$?z=il\*-&(eh-tM)\ a +£\tk-ei=i\']. , > 

yj!k(x) =e mm <Po(x) = e mm mm <po{x) 

or, in normalized form and taking into account the scaling (jl.2[) . 

(2.14) (j)>±(x) = Ce v » c 0o(x). 
The increments 

(2.15) ^:=r^=(€i-Ci=l) 

2VA " " 

formally correspond to the Yj~ from GRW and, by the aid of (|2.8|l . their distribution under P A can 
be determined if one takes into account that they have independent centered normal distributions 
with V(Z fc ) = & = £ under Q: 

P A (Zx G Ax, ...,Z n e A n ) E Q ( X{ x ieAl ,..., ZneA „ } ||^a || 2 ) 



12A3 



E)U[2v^*-(£*-^)-^|s| 2 ] 2 
E QU{X 1 6A 1 ,...,Z„6A„} / e * " 1 0o I (x)dx) 



12A51 



(2.16) 



Ai X---XA 

a 

^ ^AiX — xA„ 

([2~l"4")l and ([2~TtJ|l agree with (|23)l and (|2~TT|) . One sees that the GRW collapse process with 
deterministic times with parameters can be obtained by restricting Diosi's collapse process 

with parameter A to the appropriate discrete instants of time. 

We show now that the full GRW-process given by the parameters a and /i converges in the 
sense of finite dimensional distributions to the Diosi process with parameter A, when — > oo and 
fia/2 = A = (constant). 

Theorem 1. Let H = — + V with a bounded potential V having bounded first and second 
derivatives. Then, for all ti, . . . ,t n > 0, 

lim P2A o (0£, . . . , 0Q- 1 = P A o (0 tl , . . . , 0, J- 1 . 

Remark 2. TTie restrictions on V make sure that the Harailtonian satisfies conditions (I) and 
(II) (cf. Condition^ from the next chapter. It exemplifies that our abstract result can be applied 
to cases of physical interest. However our focus in this paper is not on greatest possible generality. 

Proof. For simplicity, we formulate the proof in one-dimensional space. The key ingredient is the 
Trotter formula (|3.8[) . which is shown in the next section. 

Similar to the observation above, we shall express the full GRW process in terms of the incre- 
ments of the Wiener process by defining 



(2.17) : = e - <(t - r »MW>* f] 

k=0 

with 

(2.18) K^t) := max{k <t}, 

j=i fi 
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setting 

~ $ 

' ll#IU' 

and performing the reweighting 

(2.19) P^(A | X) := EqOuHV^JI 2 | X) for A G ff*. 
We shall now show that this indeed yields the GRW process, i.e. 

(2.20) # = # 

where is distributed w.r.t. P^ and 4>t w.r.t. P Qj/i . In the next section we prove the product 
formula (1X8)) . Observing that ([2~T7)l is the right-hand side of (|3ll (with H = -§A + V and 
A0(:zi) = yf\xi{){x)), the proof of the theorem reduces to showing that 

Step 1: (|2~2U| holds true 

Step 2: — + V and the position operator satisfy Condition [1] of Theorem [31 
Step 3: the resulting L 2 -convergence ip^ — » ip t implies (f>^ — » 0t in the sense of finite-dimensional 
distributions. 

Step 1: Substituting ([2TT5]) into ([2T7]) as we did in (|2"7H)) . we get 

> - ^ = Ce^^W^ ^ e' f| -~ Zfc+l|2 e-^ fe+lH 



IKIU» 



fc=0 



with ^-dependent normalization factor C. Observing that, by definition (|2.15[) . 
{Z\ 6 A\, . . . ,Z n € A„} 6 5", we arrive at the following generalization of (|2.16l) 



P£(Zi e Ai, . . . , Z„ G A (1 | X) ®5P l Q (X{z ie A 1 ,...,z„ e A„}ll^„. II 2 1 *) 



03 ^ BD E Q (x We x 1 ,..^. e A.} / I II e^«-^-*l-l'e- ^^(xjda | X) 

r r n— 1 / — n 

^ / / | TT e »0-* +l -f H'e-^^^^a^da. e -«E?^ lW d ( Zl> . . . ? ^) 

JA 1 x-xA n J V ^ 

/ — n „ „ n— 1 

= \/- / / I T\e a{ - ) - Zk + 1 -^-^ z ^\-i Xk + lH cj>\ 2 {x)&x&{z 1 ,...,z n ) 

V 71" ■/AiX-xA, 7 

/ — n - - n— 1 

= \/- / / I IT e~^~ Zk+1 ^ e~» Xk+lH <f>\ 2 (x)dxd(zi, . . . , z n ). 

V K Ja 1 x-xA„ J fc = Q 

Comparing these equations with (|2.2I) . (|2.6p and (|2.5p - the ones defining GRW -, one sees that 
P^ = P Q . M so that, in particular, (|2.20j) holds true. 

Step 2: Let M. be the domain of the harmonic oscillator. According to [T3], this is invariant 
under H t and, since A s j is just multiplication by a Gaussian f (|2.12[) . p. 31) ). also under A Sj f Next, 
we choose Af := Cg° and verify the second part of Condition [T] (I)) by an explicit calculation, 
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starting with the case V = 0: 

eJ |A ^ X ^- Xx2t ~ l) <j>{x)\ 2 &x 

<3E J | ((& - 2txf - 2t) e yl ^ t - A:r2 Tk(a;)| 2 dx + 12E J \{£, t - 2tx)e VXx ^ Xx2t \ 2 \4>' (x)\ 

• :'. :<. y le^*"^ 2 * - l| 2 |<(a;)| 2 da; 

=3 / (( 2/ -2te) 2 -2t) 2 e- ij ^d 2 /|^)| 2 d a ; 



12 y -^= y (y - 2ta)V " dy|^(a;)| 2 da: 
/27rf j 

15t 2 / |0(x)| 2 da; + 12t / |0'(a;)| 2 da; + 6 /(l - e'^ 2 )^' (x)\ 2 dx 



2 



The first two summands trivially converge to for t — > 0. The last one converges to zero by 
dominated convergence because the integrand is dominated by ^x 2 \4>" (x)\ 2 . 
If V 7^ and bounded by a constant C, then by virtue of (|3.5p 

E||(A + V){A , t - 1)0|| 2 < 2E||A(A ,t - 1)^|| 2 + 2CE\\{A , t - 1)0|| 2 -> 0. 

Condition [T] (II) was verified by [18] for the case V = 0. R was again chosen to be the harmonic 
oscillator. V ^ bounded does not change the domains, so it remains to establish 

(R 2 tP,VtP) < a (U\\ 2 + \\R^\\ 2 + P) 

for some a, f3 > 0. Indeed, 

(R 2 ^,V^) = (Rij},RV^) < \\Rip\\ - WRVi/jW < i||i#|| 2 + i||i?V^|| 2 < a(||V>|| 2 + \\R^\\ 2 + /?), 

the last inequality being ensured by the boundedness of V and its derivatives. 
Step 3: We have to show that 

lim E P , ,...,%j)= E Pa (/(0 tl ,...,0 t J) 

for all bounded continuous / : L 2 (M.) n — > M. Observing that, by definition (|2.18p . there are K u (t n ) 
jumps up to time t n , which implies that, for fixed X, f(<f>t , ■ • ■ , <fit ) is ^^^j -measurable, and 
that, by (j2~l) . 

(|2~Tg|) yields 

Ep, (/(«&,...,#„)) =e q (/(^, ...,0rjii<n 2 ). 

Taking into account (I2.8[) . we are thus left to show 

(2.21) lim E Q (/(^, . . . , II 2 ) = Eq(/(&i. ■ • ■ , &JIIV>*„ II 2 )- 

Now, by triangle inequality, 

|/(C,...,C)ll<l! 2 -/(0 flI ---,0*JII^JI 2 l 

<i/(< , . . . , <ji • m^ji 2 - hvu 2 i + • ■ • , C) - » ...,0 t jniv^ii 2 

The first summand goes to in L 1 (J7) because f is bounded and ip^ — > ipt„ m L 2 according to 
Step 2. The second summand is dominated by 4C 2 ||V't„|| 2 where C is a bound for f, so, according 
to Vitali's theorem, (12.21)) is proven once we know that 

(2-22) lim |/(^,...,^J-/(0 tl ,.-.,0 t JIII^J| 2 =O 

/J.— >OG 
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in probability w.r.t. Q. For this, let e > 0. From [19], ||^|| > almost surely for every t, so, by 
cr-continuity, one can find an a such that 

(2.23) Q(lhMI,---,IKJ| >o)>i-|. 

Since \\ — > \\ijjt n \\ in L 2 , this result also allows us to find M £ N and b > such that, for all 
H > M, 

(2-24) Q(||#J,...J#J|>&)>1- J. 

Next, we choose 6 > such that, if 

\\Xi ~ 0*a II < ||Xn ~ <PtJ\ < 5, 

then 

\f(Xl, ■■■,Xn)- f{<t>tn- ■ ■ ,<f>t n )\ < e- 
Finally, again applying the L 2 -convergence iptf —> ipt, we fix N £ N such that for all /i > N 

(2.25) Q(||< -iMI,..., Il< " M < Smm(a,b)) > 1 - J. 



Applying the estimate 



< IN-yll 

~ min(||x||, ||y||) : 



Fll 112/11 

which holds true on any Hilbert space, to i = ^ and y = -0 tf , , we find that 

wr ti -<f )tl \\ 1 ...,ui-^ Ti \\<5 

outside the exceptional sets from f|2.23f) . (|2.24[) and f|2.25f) . so 

^,...,<J-/(0 tl ,...,0 t J|>e)< e . 

3. Stochastic Trotter formula 



□ 



Let H and A be possibly unbounded selfadjoint operators in the separable Hilbert space 9) with 
domains T>(H) and T>(A), respectively. In the Diosi process H = — and A is multiplication by 
t/\x. Consider the equation 

(3 x) dipt =(-iH - ^A 2 )ip t dt + Ai/ndSt 

and call its solution flow C*. The operators H and A generate the partial solution flows 

(3.2) H t := e ltH 
and 

(3.3) A s , t : = e ^-^ A -^ A2 . 

It is natural to consider the whole "anti-hermitian" part of (|3.1j) . i.e. 

(3.4) dip t = Aip t d£ t - ~A 2 -0 t dt, 

as the 'stochastic part'. Due to the damping term ~\A 2 , its solution (I3.3[) is in L 2 (f2,f)), depends 
continuously in norm on the initial data, i.e. 

(3.5) limE||(A , t -l)^|| 2 =0 

and HV'tll 2 is a martingale, in particular 

(3-6) E||^|| 2 = ||^o|| 2 . 

These useful properties can easily be verified via the functional calculus; the martingale property 
and (|3.6p also hold true for strong solutions of (|3.1I) (see [T8]). If one increased the factor i in the 
damping term, (|3.4p would still have a strong solution with \\ipt\\ 2 being a submartingale, Holevo's 
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weak solution (|12j) to the correspondingly modified equation (|3.1[) would have the same property, 
but Mora's result about strong solvability ([H]) would not apply. Without any damping term, the 
solution would in general not be in L 2 (f2, fj). Of course, if one lets A be, for example, the position 
operator, it could still be defined pointwise. 
We assume throughout this section 

Condition 1. 

(I) There exists a dense subset M C Dom(H) n Dom(A 2 ) of Sj which is invariant under H t 
and A S}t , and another dense set Af C Dom(H 2 ) n Dom(A 4 ) such that for all <j> £ N 

limE||iJ(A o , t -l)0|| 2 = 0. 
t— >o 

(II) There exists a self-adjoint positive "reference operator" R such that, with 

G:=-iH- X -A\ 

• Dom(R) C Dom(G) H Dom(G*), Dom{R 2 ) C Dom(RA) 

• {ip £ Sj | \\i(j\\ 2 + \\R^\\ 2 < 1} is compact 

• (R 2 ip,G^) + \\A^\\ 2 < a(\\ip\\ 2 + \\Ril>\\ 2 + f3) for some a, (3 > . 

(II) is taken from [T5] and guarantees that p.ip has a so-called "i?-strong solution" for any 
regular initial value. 

The following theorem is a version of the Trotter formula for the stochastic Schrodinger equation. 

Theorem 3. Let X\, X2, . . . be independent exponential random variables, which are independent 
of the Wiener process £. Define 

X 

(3.7) nJt) := max{k £ N : V — < t}, 

3=1 ' 

i.e. the corresponding Poisson process with intensity fi, and 

n v 
k=l r 

the time of the n'th jump. 

Moreover, assume that Conditional) is satisfied. Then the following formula holds true in the 
weak L 2 -sense: 

K M (t)-l 

(3.8) C ,tip= lim ff f _ r , . TT A K k+iH x k+1 ip 

U— >00 ^ V - 1 A M' M u 

fc=0 

//, in addition, (II) is satisfied and M. fl Dom(R) is dense in Sj, the convergence is strong. 

Remark 4. Formula (|3.8p should be read in the following way: A time- continuous evolution, 
given by H t with time parameter t, is interrupted at random times by jumps, described by Ak t+i . 
The jumps are Markovian, there sizes are independent and do not depend on the waiting times. In 
the proof, the indices of Ak k±i will be treated on the same level as the time parameter; however, 
this is only a technical trick. 

Proof. W.l.o.g. we first construct the Wiener process and the i.i.d. sequence on separate proba- 
bility spaces (f2^,^,Pj) and (nx,5x,Px), then form the product space 

(fi.&P) := (fi c x ftx,^®Sx,P£®Px) 

and get £ and X as the appropriate projections. 
We define the shorthand notation. 

(3.9) <7m(*) : = TT Hxj^A kh±1 . 

Lo " * " 
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Our key idea is to decompose the total error g^{t) — Cqj into the sum over the local errors over 



time intervals 



k fc+i 



Therefore we rewrite g^(t) as a telescopic sum: 

K u (t)-l 



(3.10) g„{t) = l+ \ A ^h±lHx h+1 -1 .9 M (T fc , M ). 



If we choose an initial value ip £ M, u > and /i € N, then Conditional) ensures that for 
fixed s > and variable t > s one can apply the ltd formula to As^H^^^g^u)^. We get 



A s jH c{f _ s]gil (u)i; = gu,(u)ip-i J cA s , x H c(x _ s} Hg ll (u)ipdT + J A s ^ T AH c{T _ s) g tJi (u)ipd^ T 

1 /"* 

The formula remains valid for c = — , s = k and t — k + 1. Substituting it into (|3.10[) . we get 



%(t)-i fc+i 



g^tjtp = ip-i Y L x k+iAh )3 H Xk+1 ^_^ ) Hg M (T kltl )^ds 



k=0 
K fi(*)- 1 - fc + 



(3.11) 



12 J k A ^s AH x k+l {s-^)9^{Tk tlJ )^s 



k=0 



2 51 J A ^,s A2 H Xh+l[s ^K ) g^{T k ^)^<ls 

k=0 ~ 



=:il>+I?(t)+I$(t)+Ig(t). 

We establish the weak limit of the above expressions. Note that the operators in front of the g^'s 
will converge to the operators H, A and A 2 appearing in (|3.1|) for n — > oo. 

Our first step is to show equicontinuity in a weak sense in L 2 (£l, Sj) = L 2 (fl^, C)<E)L 2 (il x , C)®$), 
i.e. we choose £ € L 2 (f^,C), F G L 2 (ilx,C) and <fi € A/", consider the complex-valued functions 

(3.12) /^(i) : =EYCv^(*M 

and show that they are equicontinuous as a sequence in /i. Then we will conclude convergence for 
H — > oo and argue that the limit actually solves (|3.1I) . 

Since H t and A s , t preserve the L 2 (f2f, f))-norm (see (|3.6[1 ), 

(3.13) ng^tW = MdM^II 2 = K x yj\\ 2 = ||vi! 2 

holds, so the terms gn(-)ip, viewed as mappings from some bounded time interval [0, T] to L 2 (fi, 
are uniformly bounded for any %p G .f). Therefore, in order to get weak convergence of gn{-)ip for 
/i — > oo, we need only consider £, Y and from dese subsets of the corresponding Hilbert spaces; 
in particular, we restrict to 

<j> e N. 

Since L 2 (il, F>) is separable (L 2 (f^) according to the Wiener chaos decomposition ([2]), L 2 (Q, X , C) 
by construction, Sj by assumption and the whole space L 2 fl,Sj due to its tensor structure), even 
an appropriate sequence 

(3.14) [YiCi<h)l&* 

suffices; for simplicity of notation, we can take normalized vectors. 
In the following, we shall repeatedly apply the inequality 

(3.15) ||xi + ■ ■ • + x n \\ 2 < n{\\ Xl \\ 2 + ■■■ + K|| 2 ), 

which holds true for any x\, . . . , x n from a vector space with norm || • || and n € N, to the norms 
on H , L 2 (fl,C) and L 2 (il,S)). Moreover, we shall often employ the integral form of Jensen's 
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inequality 
(3.16) 







2 




. dx 


< \b-a\ [ 


J a 




J a 



According to the Cauchy-Schwarz inequality and (|3.13p . every sequence is bounded uni- 
formly in t. As for the equicontinuity, we recall definition (|3.12j) and formula (13. lip and, using 
the Cauchy-Schwarz inequality, we compute 

(3.17) |/ P (t)-/ P ( S )| 2 <E|KC| 2 -lEK^^^-^^)^! 2 ^ 3^E|(0,/_f(i)-^Gs))| 2 . 

3=1 

Now we will show equicontinuity for the three summands in (I3.17[) separately. Recalling their 
definition (|3 . 1 1 1) . we start with 

E|<0,If(i)-If( S )>| 2 

Kn(t)-l -*±1 

=E|<0, J2 // X k+1 A^^H Xk+i(s _^Hg^{T k ^)^As)\ 2 

vW-i .*±i 
=E X E £ | £ ^+1/, " ( H x k+x ( s -^HA^J, g)1 (T k>ll )^ds\ 2 

k=K^(s) M 

(3.18) ®^ ExMt) _ K ^ s)l J- X 2 k+1 E e \j^ (Hx^^HA^g^T^tfdsl 3 

< -EjcMt) -*„(«)! ^ X 2 +1 / E^KF^^,)^,^,^^)^}! 2 ^ 
p , / \ J — * p 

k=K^(s) fx 

<-Ejc|K M (t)-K M (a)| £ X 2 +1 / " E s \\H Xk+iis _, } HA h ^\\ 2 \\9ATk^)M\ 2 ^ 

where we used again Cauchy-Schwarz in the last line. Recalling definitions (|3.3p and (|3.9p . we see 
that the two norm terms within each summand are independent with respect to Pg, so we can 
separate them and then use the isometry (I3.13P and the unitarity of H t : 

E\{cf>,I?(t)-I?(s))\ 2 

re^(t)-l 

\g^,{T k ^)tp\\ 2 ds 



(3.19) ' 



fc=K M (s) 



fe=« fJ (s) 

Due to definition (|3.3p and the stationary increments of the Wiener process, the integral simplifies 
to 



(3.20) 



E$\\HA k Jfds = / E\\HA , s (l)\\ 2 ds 



• ; 2/ Z.||//(J (I .> l)o!|-(ls + 2 / E||i?^)|| 2 ds< — 



where we used in the last inequality Conditional), namely that E||iJ(ylo,s — 1)0|| 2 is bounded 
since it converges to zero. Substituting into (|3.19p and using Lemma[T](a), we get equicontinuity 
oflf(t): 

m<t>ji(t) - w>i 2 < c(V\^\ + (t s f) 
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Starting again like in (|3.18[) . we proceed with 

E|(0,^(i)-^( S )>| 2 =E X E C | Yl J h " (Hx k+1 ^)^ s ^9ATk^^ s \ 2 . 

Now we interpret the sum of the integrals as one single Ito integral, to which we apply the ltd 
isometry: 

E\(4>,lZ(t)-I»( S ))\ 2 =E X ]T y/ E^\(H Xk+i{ ^^A^ s A^g^T k ^^\ 2 ds. 

Similarly to our treatment of the 1\ term, we arrive at 

E\(4>,lZ(t)-lZ(s))\ 2 < -E|K^(t) - k p («)| <C\t-s\. 
A* 

In the second inequality, we have used the fact that — /c^(s)|, the number of jumps during 

[s,t], is Poisson distributed with mean value fi\t — s\. 

The last term in (|3.17p can be estimated in the same way. One ends up with 

E\{cl>, %(t) - I^s))\ 2 <§EM<) - K ,(s)) 2 < C\t - s\. 

Thus, equicontinuity of is proven and, along the way, one gets the estimate 

(3.21) E|(0, (ff M (t) - g,(s)m 2 < C(V\t^V\ + (t - sf). 

Consequently, there exists a subsequence (/^)fcgpj of (^) m <=n such that EY"i£i (<fii , g^ (0^) converges 
uniformly in t for k —> oo (recall definition (|3 . 1 2[) and the choice (|3.14[1 ). a subsequence of (//jf.) 
such that also El^^^iO)^) converges and so on. For the diagonal sequence (^) obviously 
EYjC,j((f)j,g tl k{-)%l}) converges for all j e N uniformly in t. From the isometry (|3.13[) and the 
Cauchy-Schwarz inequality, one infers that 

L 2 (Sl,fi) 3 7]^ lim E(r],g n k(t)t/j) e R 

fc— >oo fc 

is a continuous functional. This implies, via the Riesz representation theorem, the existence of a 
process tpt with 

lim EYC{4>,g Mt (t)ij)=EYC{<f>,^ t ) 

for all Y, C, (j>. Since the limit was obtained by the Arzela-Ascoli theorem, it is uniform in t € [0, T] 
if Y, C, 4> are fixed. 

In order to show that the constructed ipt is actually the solution of (|3.1D . the convergence of 
the summands in the weak form of (|3.1ip . i.e. in 

3 

(3.22) efc(0, g^i m) =efc(0, + EY ^^ T f k (*)>> 

1=1 

to their counterparts in the weak version of (|3.1j) . i.e. in 
(3.23) 

EFc(</>, v*> = EFc(0, V) - «Fc f (H&^ds + ETc f (A4>, V^Rs - i ;T7 f {A 2 4>^ s )ds, 

Jo Jo A Jo 

has to be checked. For simplicity, we will write fj, instead of but recall that we are dealing with 
a subsequence. The left-hand side was just settled, so we start with the I = 1 term: Applying 
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Cauchy-Schwarz inequality, we arrive at 



\EY{(<f>,I?(t))-EY( / {iH^ s )ds\ 



<CE|(0,/f(t)) - / (iHcj>^ s )&s\ 



K„(t)-1 

k=0 



^A^Hx (a _ h) H gil (T fe>M )V)d S | 2 



(3.24) 



+ CE 

+ C*E 

+ C*E 

+ CE 
+ CE 



fc=0 J 71 

fc=0 ^ M 
t 



(H(f>, .g M (T LMsJjA1 )V')ds|" 



(#0, (.9m( t L^J,m) ~ fM( s ))V')ds| 2 
(i/0,. 9AI (s)V-^)d S | 2 



: Ew 



In the second step, we have substituted definintion (13. lip , expanded the difference in a telescopic 
sum and applied the inequality (|3.15p . 

Again, all the summands have to go to for /i — ¥ oo. We start by expanding the square in Jf: 



K M (t)-l 

J?(t)E x E {X J+ i-l)(X k+1 -l)E^ 

j,fc=0 



Estimating the integrals by ^ as we did in (|3.19[) and (|3.20[) and writing the sum as a square 
again, we continue 



C 



«„(*)-! 



c 



(3.25) 



Jf(«)<-2 Ex E (X k+1 -1)\ =-a(E...l Mt) < 8/lt +E...l Jt(t(t)>6/lt ) 



C" 



< — (E . . . l„ M ( t )<6pt) + (E . . . l K(1 ( t )>6/ii)'' 



Taking into account the independence of the X k , 



6/it-l 



i(E . . . l K(i(t) < 6Mt ) 2 < ^E E (X k+1 - If < -> 0, 



and using Jensen's inequality for E, (|3.15[) for the sum and the crude estimate < fit which, by 
definition (j3~32l . holds true for fc < n^{t) and, via (|535]) . implies (X k+1 ~ l) 2 < 2^ 2 i 2 + 2, 



re M (t)-l 



^(E...l Kfi(t)>6Alt ) 2 < iEK„(t) E (^fe+i-l) 2 l KM ( t )>6M<CIE KAI (i) 2 lK tI 



(t)>6/it) 



fe=0 



ON A STOCHASTIC TROTTER FORMULA 



13 



so that Lemma [T] completes the proof that (|3.25|) and thus Jf (t) goes to 0. 

All the other summands in (I3.24j) are to be treated analogously, and we shall only indicate the 
points at which new arguments are needed or the conditions enter: 

One gets 

(3.26) J£(f)<C sup E\\H(A 0tt - l)<p\\ 

o<t<J- 



2 



which goes to by condition (I), and 

(3.27) Jg(t)<CE\\(H Xk+1 -l)H<p\\ 2 . 
Since, by assumption (I), <j> 6 Dom(H 2 ), we can calculate 

(3.28) \\{H t - l)H<t>\\ = || f H s H 2 <pds\\ < t\\H 2 <t>\\ 

Jo 

and, substituting into f|3 . 2T[) . end up with 



(3.29) Jg(t) < -EX 2 

fj, 



C 

Furthermore, 



For J^, we need the equicontinuity estimate (13. 21)) : 



(3.30) ° t 



<CJ2 I " E(^/[T^~7\ + (T k ,^s) 2 ) 

1 n J — 



fc=o ' 

Now, within each integral, 

E(T fe)P - s) 2 = YT k4l + (s - f < A + J_ < 

fj, fj, 2 fj, 2 fi 

Ey/\T kill -s\ < (E(T k> „ - s) 2 ) 1 * < Cfx-i 

and, substituting back into (|3.30p . we end up with 

■TO < C^i -*-0. 

Jg (t) —> by construction of ipt as the uniform weak limit of g^(t). 
For the j = 2 summand in (|3.22[) . 

|EFc(0,/ 2 n W) -EFc^W^ s )d6l 2 < C|EC J...d^ s \ 2 . 

This can be reduced to an ordinary integral and treated like the preceding ones if one restricts 
w.l.o.g. C to the total set of iterated stochastic integrals [TT] defined via 

(3.31) M := l,M l+1 {t) := / a(s)M l (s)d( s 

Jo 

with indicator functions a and applies the ltd isometry 

E( / /d& / gd6) = e| / s d*. 

The last summand of (I3.22j) goes through like the first one. In place of ©, the terms 

E||(A 0ii -l)A^|| 2 

with j = 1 and 2 will come up, and one needs the condition <fi € Dom(A ) in order to conclude 
from p.5p that they go to 0. 
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Consequently, the constructed process ipt weakly solves equation (|3.1[) . As far as its solution is 
unique, it agrees with the one obtained by Holevo f|12|). We shall review his uniqueness proof in 
the appendix. The same method also allows to start with any subsequence gjiit)^ and extract a 
sub-subsequence converging weakly to the same ipt, so g^{t)ip — 1 ipt holds true without restriction 
to a subsequence. 

Compared to the right-hand side of (|3.8[) . the first factor is missing in the definition (|3.9p of 
g^it). Of course, both versions have the same limit because one single factor approaches the 
identity for /i — > oo: 

K M (t)-l 

EY((<t>,(H t - T ^ thit J] A± *±iffx fc+1 - 9fl (t))j>) <CTE||(g t - r ,„ m „ - 1)<P\\ 2 
fc=o " M M 



c 

where we took into account that t — T K (t),^, the waiting time for the last jump of before t, has 
exponential distribution with parameter fi. 

Since lim n _ i . 00 E||g M (i)-(/>|| 2 = \\tp\\ 2 holds true and, under the conditions from [18] . ip € AiHD(R) 
also satisfies E||V>t|| 2 = \\ip\\ 2 , m this case the convergence g^(t) — > Co,t is even strong. Since Co,t 
as well as all operators g^ have L 2 (fl, $))-norm 1, the strong convergence on the dense subset 
Ai n D(R) implies the convergence on the whole of Sj. □ 

Remark 5. It is not essential whether we group the factors in the order Ak k±± H x k+1 , as we did 

in the telescopic sum expansion (|3.10[) . or in the order H x k+1 Ak k+i . If we had decided to do the 

latter, we would have had to require 

\\A j (H t - 1M| < Ct for j e {1,2} 

instead of 

\m\m\H(A ot - 1UII 2 = 
t— >o 

in condition (I). In both cases, one can get the product formula in the form (|3.8D or with Ak *+i 
coming before H x k+1 . 

I still remains to provide the following technical lemma. 

Lemma 1. As in Theorem^ let Xi,X2, ■ ■ ■ denote independent exponentially distributed random 
variables and let K^it) be the "number of hittings" 

( k X 

(3.32) K M (t) :=maxUGN:^]— <t 



before t. Then 



a) — 



k=K^(s) 



< C(yJ\t - s\ + \t - s\ 2 ) uniformly in fx e N 



b) lim Ex K „(t)>6pt^( < ) 2 = f or ever V *• 



Proof. We start with the case s = 0, decompose 

«!•(*)! E =72 E [--- 1 «M(*)<6Mt] +72 E [--- 1 «m(*)>6^] 



(3.33) ^ 



fc=0 J " 



(you will see where the "6" comes from) and find 

6/it-l 

— E. . . l K(l(t )< 6Mt < — E E6 ^+i ^ C<2 - 
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As for the second summand in (|3 . 33[) . we know from the definition of that at least all Xk < /it, 
so 



(3.34) 



72 E - ■ • 1 K fI (t)>6^ ^ t 2 E|K. M (i)| 2 l K . M ( t )>6 M ( . 
/i 



K^(t) is Poisson distributed with parameter /it and the important point which we must not hide 
by too crude estimates is the exponential decay of P(k m (£) > 6/ii) in fit, therefore we compute 

- jU / 



fc! 



(/tf)*- 1 



^E. k pji =0,, ^(. _E 



=Ct 2 e^V lit J2 



E 



fc=([6/ii"|-l)V0 

(»•()*' 



fc! 



E 



fc=(|"6/rf"|-l)V0 



fc! 



fc! 



fc=([6/itl-2)V0 fc=(|"6//.t"|-l)V0 

Writing 

{6/tt} := ( [6/xtl - 1) VO, 
applying the Stirling formula and observing that 



(3.35) 
we obtain 



^ < I if {6/rf} > 1, 



{6/it} ~ 3 



- OO 

— E . . . l„„ (t j >6/1( <CT 2 e-^(( M t) 2 + M t) 



fc! 



< 



> 1 



(3.36) 



< 



Ct 2 e-^(( M t) 2 +M*)^?Er=o (t^j)" fOT { 6 M0 
Ct 2 ((fit) 2 + fit) for{6/tt} = 

Ct 2 ((/j,t) 2 + fit) for{6/it} = 

'ct 2 ((/rf) 2 + /if) (|) {6Mt} for {6/it} > 1 



Ct 2 ((fit) 2 + fit) 



<Ct 2 ({6 M t} 2 + {6 M t} + l) (|) 



for {6/xf} = 

\ {6/it} 



< Cf 2 . 



This proves part a) for s = 0; note that, in going from (|3.34p to the second last step in Q3.36p . we 
have also proven part b) along the way. 

In case that s ^ 0, one would be inclined to say that 



1 

7- 



E 



K M (t)-l 

Mt)-K M (*)| X k+1 

fe = K^(s) 



= ^jE 



fc=0 



because the Poisson process k m has stationary increments. In doing so, one would, however, neglect 
that the "first" X summand decomposes into 



re M (s)-l 



Kp(s) + 



i = (E ^fe+i-Ms) + (Ms- J2 X k+1 )=:Y + Z, 



k=0 



k=0 
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where ^ is the waiting time for the last jump of k m before s and ^ the waiting time for the next 
jump after s. A correct use of stationarity implies 

K„(t)-1 « M (|t-s|)-i 

(3.37) -^E\^(t) - Kfl (s)\(Z 2 + £ X 2 k+1 )^-E\ K ,(\t-s\)\ Xl+x<C\t-s\\ 

" fe=K fI (s) + l ' k=0 

Taking into account that Y and Z have independent exponential distributions and that 
x l ( a )+i < 2Y2 + 2 ^ 2 , we end up with 

—E\ Kfl (t) - Kfi (s)\ X l+i 

fe=fs„(s) 

1 . M*)-i 

<-E| KAl (i) - ^( S )|(2F 2 + Z 2 ) + -jE|« M (t) - K„(s)|(Z a + Y, X k+i) 

k=n li (s) + l 

<\(E\^(t) - ^(s)\ 2 )i(E(2Y 2 + Z 2 f)i + C\t - s\ 2 < C(^\t^J\ + \t - s\ 2 ) 
A* 

where, in the second last step, we have used Cauchy-Schwarz and (|3.37[) . □ 

Appendix: Uniqueness of the weak solution of p.l[) 

We recall Holevo's proof (p. 487 of [H]) that weak solutions of dissipative stochastic equations 
are unique under the growth condition 

(3.38) MM? < ||^ || 2 e ct for some c, 
provided that 



— iH + —A 2 is maximally accretive, 

and explain one essential step in the his proof which remained unclear to us in |12] (see below). 
For simplicity, we formulate things in terms of our equation p. II) instead of his more general 
framework. 

Suppose that ip} and ip 2 both solve p.l[) with the same initial condition i/jq and satisfy (I3.38[) . 
Defining St := ijj} — ip 2 , we have to show EMi(t)((f>, 5t) — for all 4> in a dense subset of and 
iterated integrals Mi as specified in p.31[) . The weak form p.23[) of equation p.l[) gives 

(3.39) MoW(0,5 t ) = (0,5 t )= [ (Ac)>,6 s )d( s + [ ((iH-lA 2 )<j>,8 s )ds 

Jo Jo z 

and, according to the product formula for Ito differentials, 
(3.40) 

Mtfj(<t>,6 t )= [ M,_i(a)a(s)(0, <5 s )d£ s + f Ws)(A(f>, S s )d£ s + [ M^[s)((iH — ]-A 2 )4>, 5 s )ds 
Jo Jo Jo z 

t 



Mi_i(«)o(«)<^,* a )dn 

for I > 1. Taking for granted for a moment that the stochastic integrals have mean value and 
defining mi(t) := EMi(t)6t, these equations imply 

/"* 1 

(3.41) (0, mo(t)) = J {(iH - -^ 2 )0, m (s))ds 
and for I > 1 

(3.42) (&m,(i)) = / (o(*)^,m l _i(s))d«+ / ((iff - ^A 2 )^, m z (s))ds. 

Jo Jo z 

and one can inductively prove mi = 0: Starting with I = 0, we somehow have to bring the two 
terms mo in equation p.4ip together. For this sake, we take the Laplace transform, perform 
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partial integration on the right-hand side und thus get rid of the additional integral in front of 
one mo: 

e- As (0,m o ( S ))d S = f e~ Xs f ((iH - 1 A 2 )^, m (r))drds 
Jo Jo * 

=1- J\- Xs ((iH -±A 2 )<p,m (s))d S 
Multiplying with A and dragging out the scalar products, we get 



(3.43) ((iH-^A 2 -\)<j), j\ 



mo(s)ds) = 0. 



Now let A > 0. Since — iH + \A 2 was assumed to be m-accretive, iH — \A 2 — A is surjective and 
p.43[) implies f Q e~ Xs mo(s)ds = for all A > and thus m = 0. In the induction step Z — 1 Z , 
the first integral in (|3. 421) is according to the induction hypothesis, so mi = follows from this 
equation along the same lines as mo = 0. 

We still have to prove that the sum of the stochastic integrals in p.40[) . i.e. an expression of 
the form 



J t := f Mi-i(s)X s dt a + [ Mi(s)Y s dt s 
Jo Jo 



with 



(3.44) sup EX 2 , sup EY 2 < oo 

0<s<t 0<s<t 

according to p.38p . has mean value 0. Note that the integrands are products of complex L 2 (Vt)- 
processes, so the usual L 2 -theory of stochastic integrals, providing their martingale property, does 
not apply. Holevo claims without proof that already 

sup E\Mi- 1 (s)X a \, sup E\Mi(s)Y s \ < oo 

0<s<t 0<s<t 

implies the martingale property of Jt, but we do not see the reason. 

Our argument is as follows: We compute EJ t , making use of the additional facts that Jt € i 1 (57) 
(because all the other summands in (I3.40j) are) and that the Mi are much smoother than typical 
I/ 2 -processes: 

Define the stopping times 



<r n := inf{i > : |M;_i(i)| > n or \Mi(t)\ > n}. 



Then 



Jt = f M^ 1 (s)X s dts+ f M^(s)Y a dC 
Jo Jo 



+ f (M i _ 1 (s)-Mf_" 1 (s))X s d6+ / {Mi(s)-M^{s))Y s d£ s . 
Jo Jo 

The first two integrands are majorized by n|X s | and n|V s |, so, in view of p.34p . the integrals have 
mean value zero. The last two integrands are zero before being stopped, in particular on {a n > t}, 
so the Riemann sums approximating the integrals in probability and consequently the integrals 
themselves are zero on this set. Summarizing, 

EJ t =E X{an < t} Qf (M,_i00 - M^\{s))X s dZs + J {Mi{s) - M^{s))Y s d£,)j 

=E X{an <t } Jt-^X{* n <t} [ AfJL» 1 (s)X i d^,-E X { <rB <t} / M^( S )Y s d^ s 

Jo Jo 

Now all we need is a sufficient decay of P{{r n < t}) for n — > oo, namely 
(3.45) P({a n <t}) = (\ 
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Then ~Rx{a n <t}Jt —> for n — > oo by the dominated convergence theorem and Cauchy-Schwarz 
inequality and Ito's formula would give us 

\^X{a n <t} [ M^^Xsda 2 <nWn<t}) 2 E\ f M,*» 1 (a)JT.de.| a < ^n 2 f E\X s \ 2 ds -> 
Jo Jo 71 Jo 

and the same for the last summand. According to Doob's inequality, 

P({<ro < *}) = P({ ™P 161 > n}) < and 

P(R < t}) < P({ sup |M,_i(a)| > n}) + P({ sup > n}) < E l M '-i(*)l 3 + E|M ' (t)|3 

o<s<t o<s<t n 

so that (13.451) is proven once we know E\Mi(t)\ 3 < oo for all 1. We prove this inductively, starting 
from E|Mi(i)| 3 = E|£ f | 3 < oo and using the Burkholder inequality (e.g. theorem 73 in [31]) in 
order to write E\Mi(t)\ 3 in terms of E|M;_i(t)| 3 : 

E|M,(t)| 3 =E| / a a M,_i(a)d^ a | 3 < CE[ f o,Afi_i(»)df„ / o.Af,_i(«)d^.]i 
Jo Jo Jo 

3 

(/" \ 2 (1) _ /" 

/ a 2 Mi_i(s) 2 ds ) f\ 7:; / a 3 M ; _i(s) 3 ds < C^E sup |M z _i(s)| 3 

JO / JO 0<s<t 

<8Ct*E|Afi_i(t)| 3 
[•, •] denotes the quadratic variation and can be computed according to 

[ / x s dt s , I x s &Q = I x 2 ds, 

Jo Jo Jo 

in (1) we used Jensen's inequality and in (2) again Doob's inequality. 
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